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Abstract 

In paper approach of double complex SUSY-transformations with not coincident complex energies of 
transformation is developed, allowing to deform given real potential Vi (with keeping of its shape real) 
with obtaining exact solutions. The explicit solutions of the deformation of shape of the potential, its wave 
function at arbitrary energy, not coincident with energies of transformation, wave functions at the energies of 
transformation are obtained, condition of keeping of continuity of the solutions and isospectral condition are 
determined. Using a rectangular well of finite width with infinitely high walls as the starting Vi with discrete 
energy spectrum, by the proposed approach new types of deformation of this potential with deformation of 
the energy spectrum as a whole have been obtained. The new potential contains the rectangular well as own 
partial case (with simultaneous transformation of the shape of this new potential, energy spectrum, wave 
functions of all bound states, wave function at arbitrary energy into corresponding characteristics of the 
rectangular well at needed choice of parameters). Using null potential as the starting Vi with continuous 
energy spectrum, new form of reflectionless real potential has been constructed. This potential generalizes 
well-known reflectionless potential of the type Vrct{x) = j4^(1 — 2sech^yla;), allowing: to pull down tails of the 
potential Kef in the asymptotic regions up to zero (with keeping of nonzero depth); to pull down continuously 
the depth of the hole; to displace arbitrary along axis x the hole with its passing through zero; to create and 
to increase the second hole, transforming Vrcf into double-well potential; to control continuously and simply 
the asymmetry of the shape of such reflectionless potential. Note relative simplicity of the found potential 
in a comparison with variety of the reflectionless shape invariant potentials. 
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Designations 

In paper we shall use the following designations. 

• The first potential Vi , on the basis of which we shall construct new potentials using SUSY-transformations, 
we shall name as the starting potential with number "i ". 

(k) (k 1) (k 2) 

• <pl and ifil , ipl — general and two partial solutions for wave function (WF or WFs in the plural 
case) of the potential Vk with the number k, where instead of index i at the bottom we shall use one index 
from the following: 

— n the number of the level of discrete energy spectrum of the potential Vk with number k, if this level does 
not coincide with energy of transformation £ (i. e. the energy, at which the superpotential is defined); 

— m — the number of the level of the discrete energy spectrum of the potential Vk with the number k, if this 
level coincides with the energy of transformation £; 

— w — the index at energy, if it coincides with the energy of transformation £ and does not coincide with any 
level of the energy spectrum of the starting potential Vk- 

• ifi^''^ — WF of discontinuous state, which we shall denote by stroke above. 

• E^''^ — the level with the number i (corresponding to WF iff'^) of the energy spectrum for the potential 
Vk- We shall number the levels of the discrete spectrum so: 

— the lowest level, concerned with the ground bound state, by number "1"; 

— the next levels, located higher and described excited bound states, by numbers "2" and larger. 

• ^i'^ — function of transformation {function of factorization), defined as partial solution of ip\''^ (for the 
bound or unbound or discontinuous state). 

• Sl'^^ — energy of transformation {energy of factorization), defined at potential Vk with the number k, 
coinciding with E^^\ 

• Wj;''^ — superpotential with the energy of transformation S^''^ and function of transformation ^■'^^ con- 
cerning potential Vk- 

• Operators Ai and Af we define so: 

Sometime, in text upper or lower index at such notations will be omitted. 
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1 Introduction 



At present time an essential progress has been achieved in study of properties of quantum systems. Here, a 
great importance is directed to development of methods for control (management) of spectral characteristics of 
these systems, realized by inverse problem approach [TJ[21[31[l]i smd by approach of supersymmetric quantmn 
mechanics (SUSY QM) H 1 [9l [lOl E] . 

In the approach of SUSY QM, ways of construction of new isospectral potentials (with obtaining of exact 
solutions for wave function and energy spectrum) on the basis of the given one have been studied the most 
deeply [12j . Darhoux transformations of the first order give this naturally, allowing to come from the given 
potential with discrete energy spectrum to new one with the energy spectrum without the lowest level (or with 
its little displacement for some potentials) [10l[9l[T3l[T4l[T6l[T71[T8] (see also reviews about these transformations 
|19j . about factorization method j211 |20j and historical paper [22j). Transition to the Darboux transformations 
of the second order gives a possibility to delete not only the lowest but arbitrary level in the energy spectrum, 
with obtaining explicit form for the new potential (without divergencies) and its wave functions [231 1211 IHl US] 
(see also [Ul [55]). In definition of function of factorization a transition from wave function of bound state to 
its general solution gives a possibility to select energy of factorization arbitrary. Such a way allows to analyze 
what happens with the potential if to change continuously a slope of wave function at selected coordinate at 
arbitrary selected level or to shift continuously this level in the spectrum. The constructed by such a way new 
potentials can be considered as the continuous deformations of the given potential. Transition to the Darboux 
transformations of the higher order [231129] and also development of methods of nonlinear supersymmetry [301131] 
give else richer algorithms of construction of new potentials with obtaining exact solutions of their spectral 
characteristics. Another way to construct more expansive types of transformations with obtaining new exactly 
solvable potentials can be reduced to introduction of more complicated representations of the superpotential 
(with keeping of the order of the differential transformations). Here, we note method of extended supersymmetry 
[32 ] 133 ] [34 l I35j. matrix variants of transformations of polylinear (nonlinear) supersymmetry [36] . methods of 
solution of system of coupled discrete Schrddinger equations [371 [35] > generalizations of the formalism of SUSY- 
transformations into 3 and multi- dimensional spaces }10[|ll[[39ll40] . Else one way to construct the new exactly 
solvable potentials (which can be considered as independent) consists in the development of new types of the 
differential transformations and their deformations (with keeping of the starting potential and obtaining of new 
hamiltonians of the Schrodinger type), where the standard differential operators (for example, the standard 
Darboux operators) appear as a partial case only. One can refer to the methods of nonlinear supersymmetry, 
the methods of hidden supersymmetry (hidden SUSY), developed in the approaches [HI [42l [43l [44l [45l [46] . Note 
papers [47l [48] with a possible last and full version of the method, unified these approaches. 

However, inclusion of more complicated representations for the superpotential leads to more complicated 
representations for the potentials, obtained by the approach of SUSY QM. So, use of the complex superpotential 
gives new varieties of the potentials, which turn out complex inevitably [321 [SHI [HI [51] . By use of the SUSY 
QM methods the exact solutions for new classes of potentials with PT-symmetry have been opened [S31 [SH [SS] 
[Ml [53 [511 [59] and properties of this symmetry have been studied [60l[6ll[62l[Ml[65l[Ml[66l[6Zl[68|. In research 
of new types of the potentials with real values, a difficulty in analysis of their bound states, in determination 
of divergencies which can cross out found new solution, increases essentially. How do SUSY-transformations 
with the complex superpotential or its more complicated representation expand a set of the early known real 
potentials, how difficult can bound states be calculated, how is the energy spectrum deformed? Answer on 
these questions will point to effectiveness of inclusion of the complex superpotential or its more complicated 
representations into the SUSY-transformations in search of new forms of deformations of the real potentials. 

The given paper is directed to study of these questions. We restrict ourselves by consideration of one- 
dimensional real potentials only. Using an approach of double Darboux transformations with complex super- 
potentials and not coincident energies of transformations as the SUSY-transformations, we shall analyze how 
they allow to construct new potentials on the basis of given one. The paper is organized so. In sec. [21 a variant 
of scheme of determination of partial and general solutions for the superpotential, assuming it and the function 
of factorization to be complex and the given potential — real, is presented. In sec. [3] a formalism of the double 
transformations with not coincident complex energies of transformations is developed. It turns out difficult to 
make conclusion about effectiveness of these transformations in construction of new forms of the real potentials 
in a general case. Therefore, for analysis as the starting potential we shall use rectangular well of finite width 
with infinitely high walls in sec. [4] as example of system with completely discrete energy spectrum, and null 
potential in sec. [5] as example of system with completely continuous energy spectrum. In sec. [6] conclusions are 
presented. 
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2 Superpotential as a solution of Riccati equation 



To clarify, which the most general form an unknown superpotential has, if we know only the potential, concerning 
with this superpotential in the most prevailing formalism (for example, as in [5] p. 275-277, [1^ p. 2922-2923, 
2925-2927, [TD] p. 19-20, at an energy of factorization, located not higher then the lowest level of energy 
spectrum of the given potential) , we inevitably come to a problem of a solution of Riccati equation, where as 
unknown function the superpotential is used. 

Let's consider a quantum system with one-dimensional real potential Vi(x) and discrete energy spectrum 
completely. One can find spectral characteristics (energy spectrum, wave functions) of this system, solving the 
Schrodinger equation: 

hi^i^\x) = (^-^ + v,{x)^^i^Xx)^E^J,^^gHx), (2) 

where hi is hamiltonian of the system, ipm'^ (x) is wave function (WF) of a state with a number m, Em^ is energy 
level, corresponding to the state with the number m and having real value, e[^'^ is the lowest level of the energy 
spectrum with the number 1. For the given potential Vi we introduce a superpotential W{x), defining it by 
such a condition (for example, according to 5], p. 287-289, if £ = e[^''): 

Viix)=W^ix)-'^ + 8, (3) 

where £ is constant, which can have arbitrary value and we name it as energy of transformation (or energy of 
factorization). In this paper W{x) and £ are complex. Let the potential Vi{x) be exactly solvable, i. e. one can 
write WFs for all states and the energy spectrum of the system with such a potential in the explicit analytical 
form. Let's assume, that we know a form of the potential Vi{x) and a general solution of wave function ip^^\x) 
(before imposing boundary conditions) at arbitrary energy E<-^\ Also we assume, that we do not know the 
superpotential W{x) and we shall find it. Condition defining the unknown W{x), is the Riccati equation. 
Let's find the superpotential, solving this equation. 

At first, let's find a partial solution of W{x). Fulfilling substitution of variable: 

X ^ z ^ z{x) = u{x) + i v{x), (4) 

with such condition: 

fix) = (5) 
ax 



where u{x),v{x) are real functions and f{x) is arbitrary else. Write: 

d du d dv d d d du dv 

dx dx du dx dv ^ du ^ ' ^ da; ' ^ dx 

With new variables the equation ([3]) has such a form: 



(6) 



.rl ^ c w2, ^ dW{u) , dW {u) dW {u) 

Vi{u) -£ = W\u) = W\u) - ^ ' - v^ . (7) 

dx du dv 



Let's introduce new function Wii{z): 



-uq — ivQ, at M < 0, V < 

, , 1 _ I -uq + ivo, at u < 0, V > 

~ ~JT7o' uo-ivo, atu>0, i;<0 

1*0 -f ivo, at u > 0, V >0 



(8) 



where u, v are real positive constants, zq is parameter, used for exclusion of divergence of Wh(z) at zero 2 = 0. 
We obtain: 

dWniz) d f 1 \ 1 1 



du du\ u + iv + Zq J (u + iv + Zq)^ {z + Zq)^ ' 

dWR{z) _ d f 1 \ _ i _ i 



dv dv\ u + iv + Zq J (u + iv + Zq)^ {z + ZqY 

1 



Wl{z) 



{z + zaY 
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One can write: 

where a and /3 are unknown arbitrary constants. Taking into account of this, rewrite equation ([7]) so: 
P A dW{u) 2/ _^.^.„,2 rfW^fl ^rfW^fll 



Fixing f[u) by condition: 
from ([TT|) we find: 



W^(x) = (a + j/3)/^(w)W^(u), (12) 



Vi{u)-£ = -Uj; — vx — 1-"/ — ^ ^Pf 



du dv du dv 

^—^dfWR ^^dfWR ,2dWR ,dWR 

UxVa + ip — Vooy/a + ip — ha/ — \- pj — — . 

du dv du dv 



(13) 



Now we use a and P so: 
Then from (fT3|) we obtain: 



af = Ux^/a + iP, Pf = Va:^/a+Tp. (14) 



df df df 

Vi{x) -£ = -UxWrs/oTTP — - VxWR^/a + iP — = -Wr^/o+TP ("x^ + ^'2=^) = 

, /df du df dv\ I df > (i} z(x) 

= -W^flV^TfT^f-f — + -f — ^-WR^a^iP^^-WR^Oi + ip 
V du dx dv dx / dx 



(15) 



dx^ 

Taking into account explicit form of ([8]) for Wr^ we find: 

{v,{x) - e) [z{x) + zo) = '^^^'^''It''^ ■ (16) 

One can fulfiU the condition (fT4)) . using a and /3 so: 

a + iP = l. (17) 
In this case, equation (jl6p transforms into the following: 

- ^'^dl^ + "^i(^) M + ^o)=S {z{x) + zo). (18) 

Now we have obtained the Schrodinger equation with the given real potential Vi{x), where z{x) + zo is a complex 
wave function (we denote it as (^^^''(2:)) for selected complex value £. So, let's write the partial solution W{x) 
for the superpotential: 

z{x) = Co(p''g\x) ~ Zq, 
fix) ^ d_z^^^^d4\x) 



dx ^ dx 



"^'■^'^ = -W)T7o--^;;^y (19) 



Wix) = f{x)WRizix)) 



1 #W(x) 



V^Pix) dx ' 
Co = const. 

Further in this paper, we shall name function ip^p {x) (its two partial solutions tp^^^ (x) and tp^'^^ {x)), defining 
superpotential by such a way, as function of transformation (according to logics in [521 US] , one can name it 
also as function of factorization) and denote it (^(x) (the partial solutions as (j)^^'^\x) and (/''^'^■'(x)). 

Analyzing logic of obtaining solution of Riccati equation and using the found partial solution W for the 
superpotential, we find its general solution: 

W(x) dx 

W{Cnew) = W + MV{Cnew), AW^(c„e^, = ^ ^ . , (20) 



^ J W{x) dx 



dx 
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where Cmw is new complex constant of integration. 

If we know the partial solution for wave function tp'^^\x) at any complex energy £, then (jl9p and (|20p 
determine the partial and general solutions for superpotential W{x) at complex energy of transformation £. A 
possibility of the energy of transformation to be complex introduces one new degree of freedom in definition of 
the superpotential W{x) on the basis of ([3]), keeping the given real potential Vi. Therefore, a set of all functions 
of transformation ipsi^x) is essentially larger then a set of all eigenfunctions ipn"^ (x), defined concerning the real 
potential Vi and the discrete energy spectrum En- 

In construction of new potentials-partners on the basis of the general solution W{cnew) the constant Cnew 
allows to deform their shapes. Another possibility to deform the new potential the imaginary part of £ gives. 
Therefore, Cnew and £ can be considered as independent parameters of deformation of the new potentials. 

Analysis: 

• The solutions for superpotential are obtained at arbitrariness in a choice of boundary conditions, which 
can be imposed on the function of transformation. Therefore, definition of superpotential only on the 
basis of bound states is partial one. 

• The function of transformation can be not only WF of the bound state. 

• The energy of transformation £ is arbitrary constant (it can be real, not coinciding with the ground or 
excited level Em\ or complex). 

• Solution proves a possibility to construct the continuous potential Vi with own bound states and 
without divergencies on the basis of discontinuous superpotential W{x) with divergencies. Therefore, 
investigations of complex superpotentials with divergencies have a sense. 

• The found general solution (|20p for superpotential, defined concerning one given real Vi, has two additional 
independent free parameters of deformation — imaginary parts of c„eto and £ (in a comparison with real 
superpotential for the same Vi). 



3 Double SUSY-transformations with not coincident energies of trans- 
formation: 

3.1 Deformation of wave function at arbitrary energy E/^, not coincident with the 
energies of transformation: Ek ^ £^^^ and Ek ^ E^"^^ 

Let's we know one partial solution ip'"f!''^\x) for wave function at energy E^^"^ (not obviously coincident with the 
levels En^ of the discrete energy spectrum) of the potential Vi. On its basis one can find the second partial 
solution (f l. '' (x) for wave function at this energy so: 

.P^(^)=.^^'(^)/-^t|--, (21) 
and we construct the general solution: 

X 

^i'\c,,x) = ^^^''\x)+c,-vi'''Hx)=ipi'^'\x)(l + Ck- f TTT^— (22) 

^ / (vie ' '{x')y/ 

where Cfc is arbitrary constant, which introduces one degree of freedom in deformation of the total wave function 
fk^'^ (we add index os state n to c). The value of Ck one can fix by use of boundary conditions for (^^^^ (we use 
this formalism from [33], p. 103-104; see also [35], p. 526). For arbitrary level En'^ sometimes it is convenient as 
the first solution ipn'^"^ to use WF of the bound state. Then c„ can be considered as parameter of deformation, 
varying of which, we transform this state into unbound one at level En \ deforming the shape of ipiP. 

In the first transition Vi V2 we find a general solution of WF for V2 at the arbitrary energy Ek (where 
Ek^£'^^'> and-Bfc ^£:(2)) so: 

(^f = 7Vi^)A^,(c^J^«(cfe), A^,(c.,J = A + H'(i)(c^J, iy(i)(c^J = -Ain0a)(c^J, (23) 
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where (p^}'' is the general solution of WF with parameter Ck, ipw} is the first function of transformation at energy 
f with parameter Cw-^ . In contrast with construction of formalism for real Darboux transformations, we use 
in formula (|23p a logarithm with complex argument. Using the formula for the logarithm with the complex 
argument: 

logz = log |z| + i argz, (24) 

we rewrite the superpotential in (f23|) so: 

W-'^ (c., ) =^ - ^ iH "^-^^ ("-^ ) I - * i (c» j) . (25) 

One can see that complex definition of the superpotential is possible and it expands essentially possibilities in 
search of new types of real potentials (in a comparison with using only real superpotentials). This formula can 
be useful in separation of imaginary part from real one in the new potentials. 

In inverse transition V2 v}'^°^\ we find the general solution for WF for Vi'^°^^ at the same energy so: 

^(i,def) ^ Ni')A+^ic^^)^(^\ck), <(c.J = + W('\c^,) = Ain^W(e^j, (26) 

where (j)w2 is the second function of transformation at energy with parameter c^^. 
From ([23]) and ([26]) we obtain WF at the energy Ek for the deformed v}'^°^^ : 

(^i^''^°^)(c,,c.,,c,„J^A^«A^f (£;fc-f(i))^«+A^(.^)A^f (27) 

where Aiy(c^,,c„,) = W'^^Hc^,,) - M^(i)(c^J. We find AW{c^,,c^,y. 

AW{wi,W2)^ V-^T ) , , . ■ (28) 



Taking into account: 



we obtain: 



^(1) (J^A - (J'A J' 

-(^W)" = (£(1) - V^)vL'l -(^ili)" = (£(2) _ Vi)^S, (29) 

(1)^(1) 

AW{wi,W2)^^- -7 . (30) 

(1)/ (i)V / (i)V (1) 



Using this, from (P?]) we find the final form of WF at the energy Ei 



k- 



(31) 



^^^^^^(c.,c.„c.J = N]^N^{Ek-E^^) .1^^ + A^^^ivi-^(g^^) " ^"f | 

If as the first or the second function of transformation, WF of the bound state at level f'^-* or f*^^^ is used, 
then one can normalize this function of transformation and obtain normalizing and iV^^' from this: 

n)_ r (i?,-£W)-V2, ati?,>5W, r (i?fe-5(2))-i/2^ ^iEk>£(^\ 

~ \ i|i?fc-f«ri/2^ ati?fe<f«; ^ " \ z|i?fe-£:(2)|-i/2^ ati?fe<f(2). ^^^^ 



Analysis: 

• Let's formulate a sufficient condition from (|3ip . when WF of the arbitrary state at arbitrary energy E^, 
not coincident with both energies of transformation f and £^'^\ remains continuous after deformation: 

(^) (^)) ' - (a;)) (a;) ^ 0. (33) 

This condition agrees with the requirement of inequality of Wronskian of two transformation functions to 
zero to provide the absence of singularities in the potential V2 , obtained on the basis of the given potential 
Vo by use of Darboux transformations of the second order in [21j . 
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• Fulfillment of condition ([55)1 provides continuity of WFs of all states at arbitrary levels Ek simultaneously, 
which are continuous for the starting Vi. 

• According to ((3T|l and p3|) , all bound states of the starting potential Vi remain bound after the deformation 
(with possible exclusion of states at the energies £(1) and £(2)) at fulfillment of one from 4 conditions: 

1) 0lV(xo) = O, 4'whxi)^0; 



2) 0S(xo) = O, (/.i}i(a:i) = 0; 

3) 0LV(xo) = O, cj)i^l{x,)^0; 

4) 0i)i(xo) = O, = 0. 



• In particularly, all bound states of the starting potential Vi remain bound after the deformation, if the 
energy of transformation S^^^^ or S^^^^ coincides with arbitrary level of the spectrum for Vi, and as the 
function of transformation the WF of the bound state at this level is used. 

• All bound states of the starting Vi remain bound after the deformation (with possible exclusion of states 
at energies S^"^^) at arbitrary energies f^^^ and E'^'^\ if one from conditions is fulfilled: 



1) 0LV(xo) = O, ct>^^l{xi)^Q- 

2) 0LV(xi) = O, (j)^^l{xa)^Q. 

3.2 Deformation of the potential 

The shape of the potential Vi is deformed so: 

V^'^"^'' (c^i ,0^,2) = Vi+ IWi (c„,j , ) , 



(35) 



AFi(c.„c.J = 2f£(2)-£W)A (36) 



^W2 



To be the new potential V^'^'^^'' real, we introduce the following condition: 

3(aFi(c„,,c„,,,)) =0. (37) 

Analysis: 

,(def) 



The potential V\ has 4 independent complex parameters of deformation: £1, £2, c^^, 0^2- 



• Continuity of the potential Vi , continuity of all WFs at all levels and energies (with possible exclusion of 
and f*^^-') after the deformation is defined by condition (|33p. 

3.3 Deformation of wave function at the energy of transformation £^^^ 

Let's find how WF of arbitrary state at the energy of transformation f'^^ is deformed. According to and 
(PS)) . one can find the general solution of the deformed WF so: 

^iy''''=NLl^A-A]^^Nil^{-i + W^^^^^^^^ (38) 



(2) 

We write the general solution for (f)w2 



t02 



and its derivative: 







r d 




TW2 W2 












\^dx 




































/,(!)/ 

JW2 






AD 

I^W2 



(39) 



- Ni\'^ . (40) 
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It needs to take into account, that the wave function ipw} before its deformation can differ from the used func- 
tion of transformation (j)w} at energy If such wave function coincides with the function of transformation 
(j}w} , then the general solution for ip^^^ can be found so: 

= ^i'l - = 1 . (41) 

Finding the derivative of this function: 

^if ' - (%T^)' - <^'^l - ^i\^''4^ (42) 



we obtain: 



(43) 



If wave function ip"^^ before its deformation does not coincide with the function of transformation (pwl at 
energy £^^\ then the function: 

= 4) (44) 

is not a solution of the Schrodinger equation with the new potential V2 at such energy. But in such a case, 
(2) 

non-zero function (^J,^ can be found so: 



(2) 

where N^i is new normalizing factor. Calculating the derivative of this function: 



(2),/ ^'"1 ^fcl <Pwi(Pwi P>^. 



we find: 

,(1), fi),/ ,(1),', (1) 



(1, dcf) (1) (2) /^(l) _ ^(2)\ 'P-'^'^k, ^^01 fk, , (1) . . 

Here, one can separate the deformation of ip^i^^ explicitly: 

(1, dcf) _ ^(1) ^(2) /^(l) _ ^(2)\ (1) , ^(1) ^(2) (^(1) _ ^{2)\ "P'^^n;' -^i-^' ,(1) 



,(1), (1),/ ,(1),', (1) 



3.4 Deformation of wave function at the energy of transformation 

Now we shall find how WF at the energy of transformation f is deformed. If to use the partial solution of 

(2) ('21 ('2I 

(/3it>2 at Lpyj2 — , then according to and the general solution for the deformed WF for potential 

V^""^ at such energy can be found so: 

^'^^ - %y^, li^^ = (49) 

Xq 

where, according to ([39]) . we have: 

^(2) ^ ^(r) ^.-.0^. - 0^. _ ^50) 

0L( 
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Calculating the integral: 

i^'h^) = j {^m^')) dx' = (ivw) j V ^ Ld., ^ 

xo xo yl>wi j 



(51) 



we obtain the deformed WF: 



Uj-^ \X )(pW2 V'^/ 

A2 (l)wl {x) 



('2') (2") 

Now we take into account, that the solution of tf wl can differ from the function of transformation (jrw!, . In 
such a case, one can find non-zero function ipw2 '^'^^^ so: 

/ J \ J,(2) (2)/ , ,(2),/ (2) 

(l.dcf) =^(2)^+ (2) ^^(2)f „A + p^(2)U(2) ^ pj{2)-V--'2Vn,2 +0^2 '/^t»2 .ggX 

S^1D2 ^^W2'^W2'rw2 ^*U)2 I rfx ) '^■!«2 ^^11)2 (2) ' > 

Taking into account ((39)) and (|40|) . we calculate the function (/s^^j and its derivative: 

Ai) (1),' ,(1)/ (1) 

(2) _ T^j{l) <Pwi'fw2 -(pWi 'PW2 
^W2 — U. 



.(I) 



'rW2 — -'"li'2 / \ 2 



and obtain the deformed WF (pw2 '^'^^^ ■ 



^(1) 



Al) (1),/ ,(!),/ (1) 

(1, dof) _ ^(1) ^(2) / c(l) _ C(2}\ 9W2VW2 ~PW2 ^W2 , (1) 

In the obtained solution one can separate the deformation explicitly: 

,(1) Al),/ (1) 

(1, def) ^ ^(1) ^(2) f^(2) _ (1) ^ ^(1) 7V(2) _ g{2)] ^^i'P^2 -0^1 y^2 ,(1) .ggx 

4 Example — deformation of rectangular well with finite width and 
infinitely high walls 

Let's consider, how the double complex SUSY-transformations are effective in construction of new real potentials 
on the basis of the given one. We shall analyze two cases: potential with discrete and potential wit continuous 
energy spectra. 

As the starting potential Vi with the discrete energy spectrum, we shall use the rectangular well with finite 
width and infinitely high walls: 

r 0, atO<x<a, 
+00, at a; < or x > a. 
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If to restrict ourselves for this potential only the discrete values of energy higher then the wall bottom, then we 
have a general solution for WF for arbitrary state at arbitrary energy Ek '■ 



unbound state with index k: Lp]^\x) — Asm{kx + 9), 
bound state with number n: {pn\x) = -i /— sa\{knx) 



/2 . . (58) 



and the discrete energy spectrum: 



2 

En = Ein^ = kl ^ '^n^ , A:„ = fcin = -n, n=l,2... (59) 

where A and 9 are arbitrary (complex) constants, WF of the bound state is normalized (is real, with accuracy 
to constant phase factor). Coordinates of nodes for WF of the bound state with number n equal to: 

xi^^^-, l^QA,2...n (60) 

and are located at the same distances, which are decreased with increasing of n. 

Now if to include into analysis the complex values of the energies for the potential ([57|) , then we obtain the 
general solution of the wave function concerning arbitrary selected complex energy £: 

ifi2{x) = a e"'' + 13 e-"'' = a e^'^e^'"'' + (3 e-^^e-'^"", c = x + ik, (61) 

where a and (3 are arbitrary complex constants, c is "complex wave vector", x = 5R(c) and k ~ S(c) are real and 
imaginary parts of constant c. Substituting this solution into the Schrodinger equation, we find the coefficients 
X and k: 

{-c' = £) - {k^~x^^m), -2xk^^{£)}. (62) 
If 3(5) 7^ 0, then from ((62|) we obtain: 

n{cf^^i^-n{£) + ^n{£f+-s{£f), ^{cf = ^(n{£) + ^n{£f + -s{£f). (63) 
If 3(f) = 0, then: 



X = 5ft(c) = 0, at £ = ^{£) > 0; 
k = 5(c) = 0, at £ = ^{£) < 0. 



(64) 



Further, as the method of the deformation we shall use the SUSY-approach with not coincident energies of 
transformation S^^^^ ^ f As two functions of transformation, we select the following fmictions: 

Mx) = ae=2^ +/3e-'=2^, ^ ' 

where S is arbitrary complex constant and ci, C2 are complex constants, defined by the complex functions of 
transformation E'^'^K 

4.1 Deformation of the potential 

Lets' find, how the potential Vi is deformed. According to ([36|) . the potential after deformation has a form: 

V^'^°^\cwi,Cw2) = Vi+ lWi{Cn,^, 0^2)1 

Substituting here solutions (j65p for the functions of transformation, we obtain: 

AFifr c ) - 2f£(2) _ A ae"^'^ + Pe-"^'' 
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The new potential Vi^'^^'' is complex. We shall be interesting in a question, whether the complex double 
SUSY-transformations proposed in this paper can give new real deformations of the potential Vi . The defor- 
mation AVi to be real, it is enough, that the function is real (we shall find only nonconstant (and nonzero) 
deformations AV): 

f{x) = (2) _ £(1)^ gC2^ + ^ e-«^^) (a* (4 - cl) e^2^ - /3* (c^ + c^) e-'^^^) . (68) 
Rewrite the function f{x) explicitly: 



fix) = (£:(2)-£:(i))(|a|2(c^-ct)e(==+=2)-_|/3|2(c* + ct)e-(^=+^2)^+ 
+ a* (3 (4 - 4) e('=2-c2)x _ Q,^* _^ e(c2-4)x^ _ 



(69) 



(71) 



Taking into account: 

C2+C*2=2 3?(C2), C2-C*= 2i 3(C2), (70) 

we obtain: 

fix) = _ £(1)) . (|c|2(c* - d) e23f(c2)x _ |^|2(^* ^ g^2K(c2)a.+ 

+ a*/3(c^ -ct)e-2'»(«2)x (c^ + cj) e2*^('=2)^) . 

On the basis of exponents exp ±23?(c2)a; and exp ±2i3(c2)a; one can combine functions sin (23(c2)a;), cos (23(c2)a;), 
cosh (2K(c2)a;) and sinh (25R(c2))a;. Write: 

(f (2) _ £{!)) . (^\a\^ic* - 4) e2«(^=)- - |/3|2(c* + 4) e-2*f(c2)x^ = 

= TVi cosh(2 3?(c2)a;) +7V2 sinh (2 5i(c2)ar) = e2it?(c2)x ^i + ^2 ^ g-2 5R(c.)x ^i - ^2 ^ 

(£(2) _ £(1)) . /3 i4 - 4) e-2'S*fe)x _ + g2i9(c2)x^ 



= Ml COS (2 S>(c2)x) + M2 sin (2 S>(c2)a;) = e^' f + ^ ) + e'^' ^('^^^^ (^-^ 

\ 2 2i J \ 2 2i 

where Ni, N2, Mi and M2 are new real constants. At the same exponents the coefficients must be equal: 
From here we find iVi, N2, Mi and M2, expressing them through a, /3, ci and C2: 



iVi = (£:W-£W)(|a|2(c5-cn-|/3p(c^+ct)), 

N2 = {£^'^-£^'^){\a\^4-4) + \p\^4+4)), 

Ml = {S(^'>-£W){a*Pi4-4)-ap*i4 + 4)), 

M2 = i{£^^'>-£^^^){-a*Pi4-4)-ap*i4 + 4)). 



(74) 



Taking into account only real values for Ni, N2, Mi, M2, we obtain condition of the real deformation of the 
potential: 

3^(^(2) - f (D) - 4) + 4))\ = 0, 

9f (f (2) - £W) i\a\H4 cl) + + ct))) = 0, 

^h£(^)-£i^))ia*Pi4-4)-aP*i4 + 4)))^0, 

sR((f (2) _ £:(!)) ia* Pi4-4) + a^* i4 + 4))) =0. 



Write: 

£(2) _ £(1) = _(c2)' + (ci)' = -SR(C2)' + 9(C2)' + 5R(ci)' - 9(ci)' + 2i ■ (-5R(C2)9(C2) + K(ci)9(ci)] 



(76) 
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Taking into account of this, we rewrite the first two equations of (|75p . separating real and imaginary parts 

_ £(1)) . + cj(^)2^ ^ _ (gfj(^)2 ^ ^^^^2^ _ 3(ci))} + 

(2) _ £(1)) . I (gfi(„)2 + 3(^)2) _ _ (5JJ(^)2 ^ ^^^^2^ ^g^^^^^ ^ g^^^^^^ | ^ 

3^(^(2) _ ^(D) . I (5jj(„)2 ^ cj^^^2^ (-3(c2) + 3(ci)) + {^{f3f + 3(/3)') (-^(ca) - 3(ci)) } + 
_ £(1)) . I (3ff(«)2 + cj(c,)2) _ 3ff (ci)) + (3fi(/3)2 + cj(/3)2) (gfj(c2) + gfi(ci)) | = 0. 

Add and subtract one equation concerning the other one: 

{di{af + 5(a)') • {5R(£(2) _ £(i)) (-3(c2) + 3(ci)) + _ £(i)) (gfj(c2) _ gfi(ci)) } = 0, 

(3?(/3)V3(/?)')-{5R(£(2) _£(!)) (_3(c2)-5(ci)) _f(i)) (5R(c2) + 3?(ci)) } = 0. 

And obtain the first two conditions: 

1) 3?(f (2) _ £(1)) (-3(C2) + 3(ci)) + _ £(1)) (S)ff(c2) - 3t{ci)) = 0, 

2) „ £(1)) (-3(c2) - 3(ci)) + _ f (i)) (gfi(c2) + n(ci)) - 0. 
If ^(£2 — £1) = 0, then from ([75)1 we obtain: 

3(ci) = 0, 3(c2) = 0, 9(51) = 0, 9(52) = 0. (79) 

At S(£i) — and '^{£2) = the exponents exp ±2z5(c2)a; in (|7T|) and ([7^ equal to 1. Therefore, instead of 
the third and the forth conditions of (j75p we must use the following one: 

(80) 



(77) 



(78) 



We obtain: 

3?(a)3(/3) = 3(a)3fi(/3). (81) 

One can see, that the coefficients a and (3 must be real only (with accuracy to the same complex 
factor), that the deformation of the potential AV to be real. 



5R(f (2) _ £(1)) 3(c2) - 3(f (2) „ £(1)) gfj(c2) = 0, 
Sff(£(2) _ £(1)) cj(c,)) _ 3(£(2) _ £(1)) gfj(ci) = 0. 



(83) 



If 3(^2 ^ ^1) 7^ 0, then adding and subtracting one expression concerning the other one in (|78p . we obtain: 

(82) 

Substituting the real and imaginary parts of £2 — £1 from (|76|) . we find: 

-5R(c2)' + 3(c2)' + 3fi(ci)2 - 3(ci)2] • 3(c2) - 2 • (-3fi(c2)3(c2) + 5R(ci)3(ci)) • ^{02) = 
= [3?(c2)' + 3(c2)' + 5R(ci)' - 3(ci)'] • 3(c2) - 2 3fi(ci) 5R(c2) 3(ci) = 0, 
-5R(c2)' + 3(c2)' + 5R(ci)' - 5(ci)'] • 3(ci) - 2 • (-5R(c2)3(c2) + 5R(ci)5(ci)) • 5R(ci) = 
= -3fi(c2)' + 3(c2)' - 5R(ci)' - 5(ci)'] • 3(ci) + 2 3fJ(ci) 5R(c2) ^(02) = 0. 
Multiply the first equation on 3(c2), and the second one — on 3(ci): 

3?(c2)' + 3(c2)' + 3?(ci)' - 3(ci)'] 3(c2)' - 2 3fi(ci) 5R(c2) 3(ci) 5(c2) = 0, 
-mc2f + 9(c2)' - 5R(ci)' - 3(ci)'] 3(ci)' + 2 5R(ci) 5R(c2) 3(ci) 3(c2) = 0. 
Adding one equation with the second, we obtain: 

[5R(ci)' + 5R(C2)'] (3(c2)' - 3(ci)') + [3(c2)' - 3(ci)'] (3(c2)' + 3(ci)') = 
= "5R(ci)' + 5R(c2)' + 3(c2)' + 3(ci)'l • (^(^2)' - 3(ci)') = 0, 



(84) 
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from here we find: 

3(ci) =±5(c2). (85) 
Now subtract the second equation from the first one in ([M]) : 

(5R(ci)' + 5R(C2)') (3(C2)' + 3(ci)') - 4 5R(ci) 5R(ci) 3(ci) 3(c2) = 0. (86) 
Taking into account ((85)) . we obtain: 



1) 5(ci) = +3(c2) : (5R(ci)V5R(c2)'-2 5R(ci)3?(ci) = (5R(ci)-5R(c2))' = 0) ^ (3?(ci) = +3?(c2)), 

2) 3(ci) = -5(c2) : (5R(ci)' + 5R(c2)' + 2 3fJ(ci)5R(ci) = (5R(ci)+5R(c2))' -o) ^ (5R(ci) = -3fi(c2)). 

(87) 

And we come to absolute coincidence between two energies of transformation. Therefore, we cannot use such 
solutions and the case ^{£2 — £1) 7^ 0. 

4.2 Deformation of wave function at arbitrary energy Ek, not coincident with the 
energies of transformation: Ek 7^ £^^^ and Ek 7^ S'''^^ 

Let's the general solution for WF (p^j}'' at the energy Ek (higher then the bottom of the well) before the 



deformation has a form: 



(p\}\x) — at sinfca; + 13k coskx, 



where ak and fik are arbitrary complex constants. According to ((3T|l . we find the form of this WF after the 
deformation: 



(^(^^■'^°^)(a;) = n\^^N^^ ■ I (4'^ (afc smkx + Pk coskx) + 



{kak — ciPk) cos kx — {k[3k + ciak) sin kx 
(c2 - ci) ae'^^^ - (c2 + ci)/3e-^22; 



(89) 



4.3 Bound states and the energy spectrum 

Now we consider the bound state with WF: 

</p^'''='^(x = 0) = 0, ^^^■■'"''\x = a) = 0. (90) 

At point a; = we have: 

,(i^dcf). _ w(i)/\r(2)./rfi^(i)_c(i)A ^rp(i)_P(2)V^,^«\ kak~-cif3k 



0, - ■{(^■'-^''')A + (^™-^-)(" + '^)jrF^iT|^ 

From here we find: 



,^ (gW-g(^))(a + /3) 
C2 (a — p) — ci (a + p) 



' C2{a-(3)-ci{a + l3) 



= 0. 
(91) 

(92) 



At a + /3 = we obtain: 

= 0, ^^k\x) = a/c sin /ex. (93) 

At a + /3 7^ we obtain: 

afe = • Tfc, (94) 

where 

_ ci g(i)-£:f) c2(a-/3)-ci(a + /3) 

' " /fc ^ k (f(i)-£:(2))(a + /3) • ^ ^ 

Taking into account Ek = k'^, we rewrite (pSj) so: 

_ ci C2 (g - /3) - ci (g + /3) C2 (a - /3) - ci (a + P) 

' fc ■ £(2)) (« + /?) (£(i)-£(2))(a + /3) ■ ^"^^^ 
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According to the previous section, if we consider only the real deformations of the rectangular well, thet the 
function: 



S{x) - £(2)^ 



a (c2 - ci) 6^=2^ - /3 (C2 + ci) e-'^^^ 
becomes real (at x G [0, a]). Taking into account, that it does not depend on k, we rewrite (|96| so: 

5(0) fc S'(O)' 
Now consider the deformed WF at the second point x ~ a: 



(97) 



(98) 



(£(1) 



.£:(2) 



a e 



■f3e- 



[ak sin A:a + Pk cos A:a) + 

.C2a\ (fcctfe - ci/3fc) COS fcg - {kf3k + CiOfc) sin kg 
' (C2 - ci)ae=2a _ +ci)/3e-=2a 



(99) 



From here we find: 



sm / 



ka ■ l^ak (4'^ - ^^'^) + S{a) {k(3k + c^ak)] = cosfca • ipk (4'^ " ^*'^) + ^(«) (k^k - cA) 



Taking into account ([94]) . we obtain: 



sin ka ■ 



Substituting Tk from ([98D, we find: 



Tfe (f^^) - fc2) + 5'(a) (A; + ciTk)\ = cos/ca • Ik^ - f^^^ + ^(a) {kTk - ci) 



(100) 
(101) 



sin ka ■ 



cos 



ci + £(1) - fc2] (^£(1) _ fc2^ + 5(0) 5(a) fc2 + 5(a) ci [ci + f - P 
/fca • fc • 1^(0) (fc2 _ £(1)) + S{a) (ci + f (1) - P) - 5(0) 5(a) ci 



(102) 



Now we have obtained equation of the energy spectrum for the deformed potential. At cos ka =/= one can 
write: 

k(A + Bk'^) 



where 



At 



A = -5(0) 5(a) (ci -5(0) 5(a) ci, 

B = 5(0) -5(a), 

C = (ci+£(i))£(i)+5(a)ci(ci +£(!)), 

D = -ci - 2£:(i)- 5(a) ci + 5(0) 5(a). 

A = 0, 5 = 



(103) 



(104) 



(105) 

we come to the energy spectrum for the rectangular well — before its deformation. For this, the following 
condition must be fulfilled: 



5(a) • (1 - 5(0)) • ci = 0, 



5(0) = 5(a) 



(106) 



5 Deformation of potentials with continuous energy spectra and new 
reflectionless potentials 

Now we shall analyze whether the double complex SUSY-transformations can give some-thing interesting in 
construction of new real potentials with continuous energy spectra. Let's consider the simplest case — defor- 
mation of null potential. Omitting external walls of the rectangular well (IST)) . we transform it into such null 
potential. Here, one can use formulas, obtained in the previous section. 
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As before, we use the functions of transformation in the form ([65|) . The deformed nuU potential y/^'^^-' has 
the form: 

^ ^ V J dx a {c2 - ci) e"^^ - (3 {c2 + ci) e-"^"^ ' 

This potential is complex. To be it real, we use real negative f*-^-* and £'-^\ and we use the conditions ([79|) . (I8T|) : 

5(ci) = 0, 5(c2) = 0, 5R(a)3(/3) = 3(a)5R(/3). (108) 

Lets' analyze behavior of the potential V2 in asymptotic regions. At limits x ±00 we obtain (in numerator 
and denominator we write only maximal components): 



V2{X — + ±00) 



|a|4|c2 - Ci|4e4Sf(^2)^ + |/?|4|c2 + ci|4e-45R(c,). 



= 0, (109) 

X — ^±00 



where gi and 32 are constants. So, in both asymptotic regions V2 tends to zero. Therefore, consideration of 
propagation of wave in this potential is a natural problem, and general solution for wave function for such wave 
in both asymptotic regions is linear combination of two plane waves exp(±i/cx), with opposite directions of 
propagation. 

Now we clarify whether V2 has divergencies and discontinuities. From (jl07p we see that this potential is 
continuous if denominator under differential in ()107|) is not equal to zero. So, the condition of existence of 
divergency at point x^ is the following: 

a (C2 - ci) e^^^'' - /3 (c2 + ci) e'"^'''' = 0. (110) 

We find coordinate of divergency: 

S^d^T^log — . (Ill) 



2c2 \a{c2 — ci) 
Using formula for logarithm with complex argument: 

logz = log |z| + i argz, (112) 

rewrite coordinate (jllip so: 



Xd^ log 

2C2 



/3(C2 +Ci) 



a (c2 - ci) 



i //3(c2 + ci)\ 



2c2 \q;(c2-Ci) 



Taking into account that Xd can be real only, we see that at arbitrary values of ci, C2, a, (3 the expression (jll3p 
is not correct. For example, if to use: 

^(fil^^sa^ ^i^^<o, (114) 

a(C2- ci) a (C2 - ci ) 

then the potential p07p has not divergencies and is continuous in the whole axis x. So, we have proved that 
with exception of some values of Ci, C2, a, (3, the potential V2{x) is continuous function in the whole axis x. 

The found potential (jllOp is reflectionless: if to start a plane wave from its left asymptotic region to the 
right, then it passes along the potential without any reflections, reaching the right asymptotic region. This is 
followed, for example, from Rule of construction of new reflectionless potentials on the basis of one given one, 
formulated in |69j (see p. 447; see some corrections in [70j one can make sure in this also, using analysis in 
p. 278-280 in ^5J), taking into account of continuity of both superpotentials on the whole axis. This potential 
is shown in Fig. [T] 

Let's consider a well-known reflectionless potential of such form (for example, which is presented by (33) in 
p. 280 [5]): 

Vrcfi = Al^-\l -1, (115) 

L cosh (arcfiX) J 

where Arefl and ajcfl are real constants (connected by condition: ^rofl = —ctrefi)- This potential belongs to the 
class of Poschl- Teller potentials (see papers |50) — in approaches of nonlinear supersymmetry with complex 
potentials, [411 |42l |43l [44] — in the approach of hidden nonlinear supersymmetry, [65l[64l[56] — in the approach 
with PT-symmetry, |71j — a generalization of Morse potential in two-dimensional space). One can see that in 
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Figure 1: Reflectionless potential: (a) — displacement of the hole with its passing through zero (the values are used: 
= -7, -20, -1000; f*^' = -5; a = 10 + ilOO, K(6) = 3); (b) — continuous pulling down of the hole close zero (the 
values are used: £^^'> = -20; £^^'> = -1, -3, -5; a = 10 + ilOO, »(6) = 3). 

both asymptotic regions it tends to nonzero limit (direct limit ^refl does not give zero asymptotic 

limits for this potential, with keeping nonzero hole), its depth at a: = is fixed by constant A. Our found 
reflectionless potential has also hole of finite depth, but one can pull down this hole continuously or displace it 
along axis x, passing zero; in both asymptotic regions the potential tends to null values. On the other side, if as 
the starting potential Vi to use nonzero constant potential with high A^ , then on the basis of the double SUSY- 
transformations one can obtain the reflectionless potential of the form (|107p . but with tails in the asymptotic 
regions as for the potential (|115p . So, use of negative energies of transformation with arbitrary values in the 
double SUSY-transformations allows to generalize the reflectionless potential (jllSp . 

If as the first function of transformation 0i instead of (|65|) to use the general its solution: 



then we obtain else larger possibilities in deformation of the null potential. We shall show shape of such deformed 
potential without intermediate calculations. In Fig. [2] one can see, how by shifting only the second energy of 
transformation £2 from zero to the value of the first energy of transformation £1 , one can transform smoothly 
(continuously) the symmetric one-well potential into symmetric double-well one. So, we have obtained the 
simple tools for control of asymmetry of the shape of such potentials. Note, that all these potentials satisfy 
to Rule of construction of new reflectionless potentials on the basis of the given one in |.69j with taking into 
account of continuity of both superpotentials, and therefore they are reflectionless (such potentials look like 
the reflectionless shape invariant potentials with scaling of parameters [7H [731 IZH [TBI 123 Etnd self-similar 
potentials [781 1791 iSOj (which can be symmetric in many cases) but here formulas are simpler essentially, and 
we do not find else in other papers such simple way for control of asymmetry of their shape) . 

It is interesting to analyze whether a class of the reflectionless potentials, obtained on the basis of group 
approach developed in [811 [82], can be expanded after takeing into account of the approach from sec. [3| for 
construction of the reflectionless potential (|107p . 

6 Conclusions 

In the paper the approach of double complex SUSY-transformations with not coincident complex energies of 
transformation is developed (see sec. [S]), allowing to deform the given real potential Vi (with keeping of its 
shape real) and its spectral characteristics with obtaining exact solutions. Note the following. 

• Appropriateness of such approach has been shown with obtaining of new forms of real potentials. 



01 (x) = aie=i^-H/3ie 

02 (x) = a2e=^"+/32e 



(116) 
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Figure 2: Control of asymmetry of the reflectionless potential (values are used: f'^' = -20: ai = 02 = 1, /3i = —132 = 
—3): (a) — continuous increase of the right sloping tail of one symmetric hole with creation of the second small hole 
in it = —0.00001, —2, —4, —6, —8); (b) — continuous increase of the second hole with transformation of the whole 

potential into symmetric double-well one (f '^-^ — —9, —11, —13, —16, —19). 

• Explicit solutions of the deformation of the shape of the potential, its wave function at arbitrary energy, 
not coincident with the energies of transformation, wave functions at the energies of transformation are 
obtained, condition of keeping of continuity of the solutions (i. e. without appearance of divergencies and 
discontinuities), isospectral condition with keeping of energy spectra are determined (for potentials with 
discrete energy spectra). 

• Using the rectangular well of finite width with infinitely high walls as the starting Vi with discrete energy 
spectrum, by the proposed approach new types of deformation of this potential with deformation of the 
energy spectrum as a whole have been obtained (see sec.|4l in particularly, we do not find such deformations 
in variety of deformations of the well, presented in reviews [U El [3l |4] ) . The new potential contains the 
rectangular well as own partial case (with simultaneous transformation of the shape of this new potential, 
energy spectrum, wave functions of all bound states, wave function at arbitrary energy into corresponding 
characteristics of the rectangular well at needed choice of the parameters of the deformation) . 

• Using the null potential as the starting Vi with continuous energy spectrum, new form of the reflectionless 
real potential has been constructed (see sec. [5]). This potential generalizes the well-known reflectionless 
potential of the type Vrcf{x) — A^{1 — 2sech Ax), allowing: 

— to pull down tails of the potential Vrai in the asymptotic regions up to zero exactly (with keeping of nonzero 
depth) ; 

— to pull down continuously the finite depth of the hole; 

— to displace arbitrary along axis x the hole with its passing through zero; 

— to create and to increase the second hole, transforming Vrd into double-well potential; 

— to control continuously and simply the asymmetry of the shape of such reflectionless potential. 

Note relative simplicity of the found potential in a comparison with variety of the reflectionless shape 
invariant potentials (with scaling of parameters; for example, see [5]). 
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